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Abstract. 

We describe the computer codes, developed at Vilnius University, for the calcu¬ 
lation of electron-impact excitation cross sections, collision strengths, and excitation 
rates in the plane-wave Born approximation. These codes utilize the multireference 
atomic wavefunctions which are also adopted to calculate radiative transition param¬ 
eters of complex many-electron ions. This leads to consistent data sets suitable in 
plasma modelling codes. Two versions of electron scattering codes are considered in 
the present work, both of them employing configuration interaction method for in¬ 
clusion of correlation effects and Breit-Pauli approximation to account for relativistic 
effects. These versions differ only by one-electron radial orbitals, where the first one 
employs the non-relativistic numerical radial orbitals, while another version uses the 
quasirelativistic radial orbitals. The accuracy of produced results is assessed by com¬ 
paring radiative transition and electron-impact excitation data for neutral hydrogen, 
helium and lithium atoms as well as highly-charged tungsten ions with theoretical and 
experimental data available from other sources. 

Keywords: electron impact, excitation, many-electron ions 


PACS: 31.15.ag, 34.80.Dp, 95.30.Ky 


Calculation of electron-impact excitation parameters 2 

1. Introduction 

In modelling both high temperature plasma (stellar atmosphere, nuclear fusion) and low 
temperature plasma, such as planetary nebulae, working material of spectroscopic and 
medical devices, one needs data on free-electron interaction with atoms and ions. For 
the consistency of plasma models, it is highly desirable that such data are calculated 
within the same approximation, using the identical multireference atomic wavefunctions, 
applying the same methods to include relativistic and correlation corrections as it has 
been done in production of spectroscopic data, such as energy levels, oscillator strengths 
and radiative transition probabilities. 

Over many years, original methods and computer codes designated to calculate 
various spectral parameters of atoms and ions have been developed in the Department 
of Atomic Theory, Institute of Theoretical Physics and Astronomy, Vilnius University 
00. Currently these codes were supplemented with our new codes computing electron- 
impact excitation parameters for ions in the plane-wave Born approximation. The 
main purpose of this development is to establish a consistent and complete set of data, 
necessary for plasma spectra modelling. Such a set will be suitable for our newly- 
developed database ADAMANT (Applicable DAta of Many-electron Atom eNergies 
and Transitions), where the main requirement is to produce data sets within the same 
atomic wavefunctions base, hence, simplifying an application of these data in modelling 
codes. In Section [2] of the present work, for the first time we describe the implemented 
calculation methods. Further, in Section |3l we describe the algorithms implemented in 
our computer codes. 

Developed methods and computer codes are equiped to calculate many-electron 
atoms and ions with open s-, p-, d- and f- shells when consistent inclusion of the 
correlation effects is necessary. It is important to perform calculation of the electron- 
impact excitation parameters using an extensive conhguration interaction basis in order 
to match them to other spectroscopic parameters, such as energy levels, transition 
probabilities, oscillator strengths) determined in the same approximation and similar 
accuracy. These codes will be employed in cases when the adaption of other more 
accurate theoretical mathods, such as R-matrix approximation (RM) or converged close- 
couling approximation (CCC)), is very difficult or even impossible due to complex atomic 
structure (e.g. for heavy multicharged ions). The plane-wave Born approximation to 
calculate the electron-impact excitation parameters was chosen as a suitable one be the 
developers of Atomic Data and Analysis System (ADAS), see [3]. ADAS needs the data 
for complex tungsten ions which can not be determined using other approximations. 

In the current work, in order to benchmark the adopted approximations and 
developed codes, we present the investigation of the electron-impact excitation cross 
sections for light atoms, namely H, He, and Li. There is substantial amount of atomic 
structure and electron-atom interaction data for these atoms, both theoretical and 
experimental ones. Most of them can be found in the NIST Atomic Spectra Database 
|1], the NIST Electron-Impact Excitation Cross Sections Database 0, and the CCC 
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database |S]. We choose H, He, and Li atoms for comparison since these are the only 
atoms having electron-impact excitation data in the NIST database [5] which usually 
contains only very reliable parameters. We have performed our calculation by employing 
both the non-relativistic radial orbitals and the quasirelativistic ones. Determined 
results are compared in Section 01 Conclusions are presented in the hnal Section. 

2. Description of adopted method 

The total cross section a, describing the interaction of an incident electron having the 
energy e with an atom, for the excitation from the energy level KqXqJq to the level 
KiXiJi is expressed as a sum of the excitation cross sections with different ranks k: 

(^{KqXqJq, KiXiJi, e) = aK,{Kf)XoJo, KiXiJi,e). (1) 


Here K denotes electronic conhguration, A denotes level number, J stands for the total 
angular momentum and e is the energy of incident electron. Contrary to the case of 
radiative transitions, cJk does not contain hue-structure parameter a, and the summation 
in (JT]) must be performed over all possible ranks k. This summation must be performed 
over even or odd n values, depending on the parity of the initial and hnal levels, including 
K = 0, if the excitation process does not change their parity. The rank n must satisfy 
the triangular condition with the even perimeter ( Jq, /«, Ji). 

The excitation cross section of any rank k is expressed using the matrix element 
iSre of electron-impact excitation operator: 


ko+k 



( 2 ) 


kQ—k 


Here oq = 0.5291772 x 10“^° m is atomic length unit. We must underline that all our 
calculations are performed in the system of atomic units. Integration is performed over 
the diherence of the electron momenta; ko is the momentum of an incident electron, and 
k is the momentum of an outgoing electron. 

Our developed methods and computer codes determine cross sections for 
transitions between diherent energy levels by employing multiconhgurational multi-term 
wavefunctions in J-coupling: 


'^{KoXoMx) = a{KoXoJo,KTLSJo)^{KTLSJo\x). 


( 3 ) 


KTLS 


Here TLS denote the intermediate and hnal orbital and spin-angular momenta for the 
pure LS'J-coupling wavefunction, and x denotes coordinates of all electrons of the 
wavefunction T. In this case, the matrix element of the electron-impact excitation 
operator is expressed by relation 
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S^iKoXoJo, q)= <^oXoJo, KTLSJo) 

KTLS,K'T'L' 

X S^{KTLSJo, K'T'L'SJi, q)a{K'T'L'SJ^, KiXiJi). (4) 


In this approximation, the excitation process does not change the term multiplicity 
S. The matrix elements for the excitation process in pure L5'J-couphng can be described 
by a product 

S^iKTLSJo, K'T'L'SJi, q) = Q^iKTLS, K'T'L'S)W^{Jo, Ji)R.{n-A, ndu q) ■ (5) 


Here Q^iKTLS, K'T'L'S) is an angular integral for the matrix element of the 
electric multipole transition operator of the rank k, which can be determined using code 
developed for the calculation of the radiative transition matrix elements. Therefore, 
the rank k must satisfy the triangular condition with even perimeter {Lq, k, Li). 
The expression for also includes sub-matrix elements of the spherical function, 
consequently, an additional triangular condition of with the even parameter {£i,K,£i) 
arises, where C and £f are the orbital momenta of excited electron before and after 
collision. The coefficient describing the bounding of the total orbital momentum L and 
the total spin S into the total angular momentum J is the same one as in calculation 
of the radiative transition matrix element: 




- 1 ) 


L'+k-J-S 


L L' K ] 

J' J ^ J ’ 


( 6 ) 


where the element in curly brackets denotes the 6j-coefficient. Rn{ni£i, nf£f, q) is a radial 
integral: 

OO 

R^{ni£i,nf£f,q) = j P{ni£i\r)P{ni£f\r)[j^{qr) - d^,o]dr , (7) 

0 


where P{n£\r) denotes one-electron radial orbital of the initial and the final state of an 
electron, 5^0 is a Kronecker delta function, and is a spherical Bessel function |8]: 


jn^qr) 



JK+i/2{qr) . 


( 8 ) 


Finally, we can define an expression for the total electron-impact excitation cross 
section: 


(^{KqXqJq, KiXiJi^e) — 


e 


(2Ji 


fco+fc 

^ J Tyrr T o Ty, 


- fco-fc KTLS,K'T'L' 


X Q^iKTLS, K'rL'S)WM: R) 


a{KoXoJo,KTLSJo) 
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X j P{n-A\r)P{niC\r)[i^{qr) - dr 
0 


2 


X a{K'T'L'SJi,KiXiLS,) 



(9) 


From ([9]) one can see that the expression for electron-impact excitation cross section 
has a similar structure and analogous angular integrals compared to calculation of the 
radiative transition probabilities. Nevertheless, the radial part in the electron-impact 
excitation cross section expression is much more complex. There is a double integral and 
sum over all possible ranks k instead of one single integral. Moreover, calculations must 
be performed for various energies of the incident electron e. These energies are usually 
expressed using the excitation threshold energies, which are determined as a difference 
between the energies of the hnal and initial levels. Consequently, there is a different 
set of free-electron energies for each excitation channel. Therefore, a special attention 
should be paid in order to simplify and speed-up calculation of the radial integrals for 
the excitation cross sections. 


3. Calculation algorithm 

At hrst, we determine all necessary angular integrals Qf^{KTLS, K'T'L'S) for the matrix 
element of transition operator. In order to avoid repetitive calculations of the radial 
integrals with close values of the parameter q, we exploit interpolation of radial integrals. 
For this purpose, we determine the minimum and maximum values of q from the initial 
information about the investigated excitation. Within these limits, a grid of 1001 q- 
values is generated. Since the functions in calculated integral are more sensitive at the 
lower parameter values, this grid of g-values is formed using a logarithmic step. Next, 
a two-dimensional grid of q and the radial variable r products is generated. A grid for 
the Bessel function (0) is determined according to the grid of the parameter q ■ r for all 
necessary rank k values, starting with k = 0. After the Bessel function grid is generated, 
all the possible integrals ([7]) are determined for given q values. These computations do 
not require too much time even if the quantities of the parameter q and radial variable 
r values are large. 

After the grid is determined, electron-impact excitation cross sections can be 
produced. In most cases, the incident electron energies e are introduced in the excitation 
threshold units. Nevertheless, there is a possibility in our code to dehne a different 
energy grid. A grid of the parameter q is determined for each incident electron energy 
in order to perform integration in ([2]). As in previous case, this grid has a logarithmic 
step. We have performed a computational experiment and have determined, that a 
desirable accuracy can be achieved if the Simpson’s rule |9] is applied for the integral 
when a number of points is as low as 25. The integrals q), which are 

necessary for calculation of the electron-impact excitation matrix elements in ([5]), are 
determined by interpolating their values according to the previously generated values. 
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Only the nearest fonr points are reqnired for this interpolation. Integration over the 
parameter q is performed for the determined matrix element. Further, the calculated 
Ki\iJi,e) values are summed over all possible ranks n as in ([T]), and the 
total electron-impact excitation cross section is determined. We are performing parallel 
computing for the cross sections of different ranks k in order to accelerate production 
of results. 

The above described calculation method is realized in two computer codes. One of 
these codes exploits the non-relativistic Hartree-Fock (HF) radial orbitals. Another one 
employs the numerical quasirelativistic (QR) radial orbitals. The performance of our 
new codes was tested on complex heavy ions, such as Hf- HF"*", Ta"*" - Ta^"*", 
and Re^"*" - Re®"*". Such complex systems require to deal correctly with both relativistic 
corrections and correlation effects. Unfortunately, there are no reliable experimental or 
theoretical data for such complex systems published so far. Furthermore, in order to 
benchmark these newly developed codes, we have performed calculation of the electron- 
impact excitation cross sections for three light neutral atoms, namely, hydrogen, helium 
and lithium. Some data for the ion were also calculated and compared with 

calculations performed using more sophisticated R-matrix method in the next section. 
Finally, we demonstrate abilities of our computer codes by presenting the electron- 
impact excitation cross sections for several low-ionization stages of tungsten ions. 


4. Results and discussions 


The NIST database [5] contains the total electron-impact cross sections for the 
transitions from the ground configuration with one level to an excited configuration 
rather than to separate levels. Our computer codes are designed to determine cross 
sections for the transitions between the individual levels of many-electron atoms and 
ions ([9]). Therefore the results determined in the present work were summed over all 
levels of the final configuration. In general case, a cross section averaging is performed 
by applying expression: 


. X^AoiJoi.Ai, Jh + l)c’‘(-RoAoi</oi,Jlj, e) 


( 10 ) 


The most of the data for comparison with our results are taken from |T0]. The 
electron-impact excitation cross sections in that work are scaled (RU-scaling) according 
to a method developed for neutral atoms in [TT] : 

(^be{Ko, Ki,e) = . n . Ki,e). ( 11 ) 

£ JD Jh 


Here B is the ionization energy of the initial level, E is the transition energy. Such 
a scaling signihcantly alters cross section values, particularly those at the energies close 
to the excitation threshold. This leads to a substantially better agreement between 
theoretical results and experimental data. Since the scaled data for cross sections are 
presented in [10], we utilize the RU-scaling for all our results in the current work. 
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Likewise in HD], we utilize the B values from [S] , and these values agree favorably with 
our theoretical results. For the transition energy E, we have utilized our ab initio 
energy values. The electron-impact excitation cross sections in [10] were computed 
using the one-configuration Dirac-Fock approximation. Therefore, in order to improve 
the accuracy of calculated data, authors introduced an additional /-scaling: 


o'BEjiKo, Ki, e) 





o’BEiKo, Ki,e) , 


( 12 ) 


where the ratio of oscillator strengths was applied. Here /sc is a theoretical (one- 
conhguration approximation) oscillator strength, and /ac is a high-accuracy oscillator 
strength from [1]. Since the multiconhguration approximation is employed in the present 
work, the /-scaling is not performed for production of our results. 

For comparison with the results from [10], we introduce a mean-square deviation 
MSD dehned as: 


MSD = 




1 ^ 

N ^ 


n=l ^ 


(^BE,f{Ko, Ki^Sn) — <Jbe{Kq^ Ki, Sn) 
<^bej{Ko, Ki,en) 


X 100%, 


(13) 


where (Tbej^Kq, Ki^Sn) stands for the data from [TD], asE^Ko, Ki, e„) is our calculated 
cross section results. Summations were performed for all electron energies presented 
in the NIST database. 


4-1. Excitation of hydrogen atom 

We have investigated the electron-impact excitation of the Is electron to the 2p, 
3p, and 4p states. As one can expect, an agreement with data from [I0| is excellent 
for all three transitions within the complete energy range. Both the HF and the QR 
results display MSD = 0.22% for these three transitions, and our results are slightly 
lower than the data from HD]. The deviations for the ls-2p excitation are larger at 
the beginning of the investigated energy range, where the deviations are approximately 
0.43% at e = 11 eV, and at the end of the presented energy range, where the deviations 
are roughly 0.82% at e = 3000 eV. In the middle of the energy range, the deviations do 
not exceed 0.12%. A completely similar situation is for the excitation to the 3p state. 
For excitation to the 4p state, the deviations at the low and at the high energy end 
increase, but this increase does not change MSD noticeably. 

In tabled] we present the electron-impact excitation cross sections of hydrogen atom 
for the cases, where it is possible to make comparison with the data produced by several 
authors. For our data, we present only the HF results because, as one can expect, the 
difference between our HF and QR results appears only in the hfth or sixth signihcant 
digit. We present the theoretical results, determined using the convergent close-coupling 
method (CCC) [6] and the experimental data from [12]. It is not unusual that our data 
agree better with the data from [10] rather than with those from [6]. Their agreement 
with the experimental data from na is within error limits. 
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Table 1. The electron-impact excitation cross sections a (in 10 cm^) for H atom 
at various incident electron energies e. 


e 

(eV) 

HF 

BEJ 

BH] 

CCC 

0 

Exp 

m 

14 

3.506 

ls-2p 

3.512 

3.849 


15 

3.920 

3.926 

4.021 

4.8 

20 

5.272 

5.278 

5.163 

5.7 

30 

6.312 

6.319 

6.187 

6.4 

40 

6.524 

6.532 


6.7 

45 

6.506 

6.513 

6.499 


70 

6.025 

6.032 

6.142 


100 

5.330 

5.337 

5.485 


150 

4.417 

4.423 

4.582 


200 

3.772 

3.777 

3.912 


500 

2.076 

2.081 

2.126 


1000 

1.242 

1.246 

1.261 


14 

0.437 

ls-3p 

0.438 

0.564 


15 

0.536 

0.537 

0.587 


20 

0.827 

0.828 

0.752 


30 

1.039 

1.040 

0.960 


45 

1.086 

1.087 

1.024 


70 

1.010 

1.012 

1.014 


100 

0.895 

0.896 

0.899 


150 

0.741 

0.742 

0.761 


200 

0.632 

0.633 

0.652 


500 

0.346 

0.347 

0.355 


14 

0.127 

ls-4p 

0.127 

0.203 


15 

0.169 

0.170 

0.209 


20 

0.284 

0.284 

0.242 


30 

0.364 

0.364 

0.325 


45 

0.382 

0.383 

0.352 


70 

0.357 

0.357 

0.355 


100 

0.316 

0.316 

0.315 


150 

0.262 

0.262 

0.268 


200 

0.223 

0.223 

0.230 


500 

0.122 

0.122 

0.125 



4 . 2 . Excitation of helium atom 

We determined the cross sections of electron-impact excitation from the Is shell to the 
2p and 3p shells. Our HF results were obtained in a following way. First of all, the 
Hartree-Fock equations were solved for the ls2s conhguration of helium atom. At next 
step, the equations for the 2p, 3s, and 3p electrons were solved in a potential of the 
frozen Is electron. The determined radial orbital basis was complemented with the 
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Table 2. Energy levels E and radiative transition probabilities A in He. 


Level 


E (cm 



A (s-i) 


m 

HF+Cl 

QR-tCl 

m 

HE-hCl 

QR-tCl 

ls2p 

169087 

169031 

169106 

1.764E-h2 

1.586E+2 

1.714E+2 

ls2p iPi 

171135 

171078 

171140 

1.799E-H9 

1.794E+9 

1.826E-h9 

ls3p 

185565 

185513 

185573 


3.093E+1 

3.312E-H1 

ls3p 

186209 

186160 

186255 

5.663E-h8 

5.604E-H8 

5.637E-h8 


transformed radial orbitals (TRO) |2llT3], which described virtual electron excitations 
and had the principal quantum number values 4 < n < 9 and all possible values of the 
orbital momentum i. Consistent with our other calculations, we employed the same 
radial orbital basis both for the even-parity and odd-parity configurations. Therefore, 
we can avoid problems related to the non-orthogonality of radial orbitals. 

The correlation effects are included within the configuration interaction (Cl) 
approach. We selected admixed configurations which had averaged weights, calculated 
in the second order of perturbation theory na, larger than 10 in the wavefunction 
expansions of the adjusted configurations ls2s, ls2p and ls3p by applying the method 
described in m ng. We can adopt such a minute selection criteria, because 
our investigated configurations have only two electrons. Therefore, the constructed 
Hamiltonian matrices are relatively small. We apply Breit-Pauli approximation to 
include relativistic effects. Furthermore, the parameters of radiative transitions and 
electron-impact excitations are calculated after the eigenvalues and the eigenfunctions 
are determined. Our QR-I-CI approximation results are obtained in very similar way, but 
this time the quasirelativistic equations [MiiniiiH] are solved instead of Hartree-Fock 
equations in order to determine the one-electron radial orbitals. 

The investigated energy levels and their radiative dipole emission transition 
probabilities to the ground level are given in table [2l We present only those levels which 
can be excited by electron impact from the ground level via the dipole transitions. For 
the He atom, our calculated results are compared with those from the NIST database. 
Those data were determined using high-accuracy calculations for the helium isotopes 
ng. It is evident from table |2] that our theoretical level energies agree very well with 
NIST data. The deviations reach only few one-tenths of cm“^. Although the relativistic 
effects for the helium atom are small, our QR-I-CI calculations give more accurate energy 
level values. Similar high accuracy is evident for the emission transition probabilities, 
when they are compared with the data from compilation [20]. The accuracy of the 
QR-I-CI results is fine even for the intercombination transition ls2p ^Pi— Is^ ^Sq. 

Our cross sections agree well with the data from [T0| for the ls-2p excitation. Some 
larger deviations, reaching few percent, are noticeable only for the low incident electron 
energies e. Otherwise, these deviations are significantly lower than 1%. Consequently, 
the mean-square deviations from all data contained in the NIST database [10] are 1.2% 
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Table 3. The ls-2p electron-impact excitation cross sections a (in 10 cm^) for He 
atom at various incident electron energies e. 


e 

(eV) 

HF-fCI 

QR-tCI 

BEJ 

dni 

CCC 

m m 

m 

Exp 

m 

m 

RV 

m 

25 

2.851 

2.733 

2.933 

2.23 



1.165 


1.2 

30 

4.566 

4.440 

4.669 

3.89 


3.75 

3.408 


3.7 

35 

5.770 

5.664 

5.875 



5.43 




50 

7.806 

7.775 

7.903 

8.19 


8.20 

8.370 


8.5 

60 

8.452 

8.455 

8.542 



9.52 

9.240 


9.5 

70 

8.799 

8.827 

8.883 



9.77 

9.660 


10.0 

80 

8.965 

9.009 

9.043 

10.60 


10.15 

9.810 


10.2 

90 

9.014 

9.071 

9.088 



10.15 

9.820 


10.2 

100 

8.989 

9.055 

9.060 

10.87 


10.10 

9.740 


10.1 

120 

8.812 

8.889 

8.876 



9.63 




150 

8.406 

8.490 

8.462 



9.18 

8.780 


9.2 

180 

7.961 

8.047 

8.012 



8.81 




200 

7.670 

7.756 

7.718 

9.05 

8.97 

8.30 

7.747 

9.2 

8.1 

250 

7.004 

7.087 

7.046 



7.63 




300 

6.433 

6.512 

6.471 



6.95 




350 

5.947 

6.022 

5.982 



6.41 




400 

5.532 

5.603 

5.564 



5.94 




500 

4.862 

4.926 

4.890 

5.54 

5.54 

5.07 

4.579 

5.2 


900 

3.330 

3.376 

3.350 

3.70 






1000 

3.096 

3.139 

3.116 



3.14 

2.903 



1500 

2.316 

2.349 

2.333 



2.34 




2000 

1.869 

1.896 

1.885 



1.85 

1.747 




for the HF-I-CI approximation, and MSD = 2.6% for the QR-I-CI approximation. The 
larger MSD value for the quasirelativistic case is caused by more than two-times larger 
deviations, compared to the HF-I-CI approximation, at low electron energies (e < 50 eV). 
The deviations from the compiled data |T0] both for the QR-I-CI and for the HF-I-CI 
results are of similar size at higher electron energies. 

We compare our results with the data from other authors in table [3l As in table H] 
we display only those electron energies, for which there are data obtained using different 
methods. Such are the theoretical results from PEI], the experimental data from 

[22] 1^ . and the recommended values (RV) from |2l]. One can see from table [3] 
that our results are closer to those from [10] than to other theoretical data. This can 
be explained by the fact that we adopt the same scaling procedure (HD as in [To] . 
Unfortunately, we can not achieve close agreement to the recommended values (RV) 
PI, especially at low electron energies. 

For the ls-3p excitation, our results agree with those from [10] slightly worse 
compared to the ls-2p excitation results. Here MSD = 3.7% in the HF approach, 
and MSD = 4.2% in the quasirelativistic approach. It is interesting that agreement for 
this transition is better at lower energies, but it exceeds 4% and 5%, correspondingly. 
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Table 4. The ls-3p electron-impact excitation cross sections (in 10 cm^) for He 
atom at various incident electron energies e. 

e HF-f Cl QR-f Cl BE, f CCC Exp RV 


(eV) 



m 

m 

m 

m 

25 

0.509 

0.494 

0.498 

0.40 



30 

1.028 

0.992 

1.002 

0.75 

0.74 

0.582 

35 

1.366 

1.313 

1.329 


1.11 


40 

1.611 

1.548 

1.565 

1.40 

1.38 

1.400 

45 

1.793 

1.726 

1.740 


1.66 


50 

1.930 

1.863 

1.871 

1.88 

1.84 

2.166 

60 

2.112 

2.052 

2.043 


2.16 

2.098 

70 

2.214 

2.165 

2.138 


2.34 

2.216 

80 

2.266 

2.228 

2.185 

2.46 

2.43 

2.265 

90 

2.287 

2.259 

2.202 


2.45 

2.274 

100 

2.288 

2.267 

2.199 

2.59 

2.44 

2.258 

120 

2.252 

2.245 

2.161 


2.41 


150 

2.156 

2.162 

2.066 


2.26 

2.041 

180 

2.046 

2.060 

1.959 


2.20 


200 

1.973 

1.990 

1.888 

2.23 

2.08 

1.805 

250 

1.804 

1.825 

1.726 


1.88 


300 

1.658 

1.680 

1.587 


1.77 


350 

1.532 

1.554 

1.468 


1.63 


400 

1.425 

1.447 

1.365 


1.51 


500 

1.251 

1.272 

1.201 

1.37 

1.29 

1.086 

900 

0.854 

0.868 

0.823 

0.92 



1000 

0.793 

0.807 

0.766 


0.79 

0.691 

1500 

0.592 

0.602 

0.573 


0.59 


2000 

0.476 

0.484 

0.463 


0.47 

0.417 


at high electron energies. Our HF-I-CI cross section values are slightly larger than 
those from [10], whereas the QR-I-CI data are slightly smaller at low energies, but 
they become larger at higher electron energies. The electron-impact excitation cross 
sections are tabulated in table 01 where we present the theoretical data (CCC) from 
[6] , the experimental data from |22l |23] , and the recommended values from [25] . The 
fact that our calculated values are larger than those from [10] brings them closer to the 
recommended data from [25] . 

4-3. Excitation of lithium atom 

For the lithium atom, we investigated the excitation of the 2s electron to the 2p and 3p 
shells. Calculations were performed in the following way. First of all, the Hartree-Fock 
equations were solved for the ls^2s configuration. Afterward, the equations for the 2p, 
3s, 3p, and 3d electrons were solved in a frozen-core potential. This basis of radial 
orbitals was complemented with the TRO describing virtual electron excitations for the 
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Table 5. Energy levels E and radiative transition probabilities A of lithium atom. 


Level 


E (cm 



A (s-Q 


NISTg] 

HF-PCI 

QR-tCI 

NIST HO] 

HF-PCI 

QR-tCI 

2 p ^Pl/2 

14904 

14923 

14900 

3.689E-P7 

3.719E-P7 

3.694E-P7 

2p ^^3/2 

14904 

14923 

14900 

3.689E-P7 

3.719E-P7 

3.694E-P7 

3p ^-Pl/2 

30925 

30924 

30915 

1.002E-P6 

0.943E-P6 

0.945E-P6 

3P ^^3/2 

30925 

30924 

30915 

1.002E-P6 

0.943E-P6 

0.945E-P6 


principal quantum number A <n < 11 and for all possible values of the orbital quantum 
number 1 . 

We selected all admixed configurations which have their averaged weights [12] larger 
than 10“^® in the adjusted conhgurations ls^2s, ls^2p, and ls^3p. Calculations in the 
quasirelativistic approximation were performed in a completely analogous way. Such an 
extension of the TRO basis and the reduction of the conhguration selection parameter 
down to 10“^®, compared to the calculation of He atoms described in Sect l4.21 was 
necessary in order to ensure the convergence of the 3p-2s transition probability. The 
determined energy levels and the transition probabilities A are presented in table [5l 
One can see that agreement of the energy level values with data from the NIST 
database is rather good. Agreement of the 2p-2s transition probability values with 
the compilation data from [20] is hne, too. For the 3p-2s transition, we cannot achieve 
such an agreement, and the discrepancy of our data is approximately 6%. 

The cross sections for the electron-impact excitation of the 2s-2p transition are 
presented in table [6l One can see that our results completely agree with the data 
from [To]. For the most of calculated energies, the deviations from the NIST database 
values are less than 0.5%. For the HF-I-CI calculations, MSD = 0.55%, and for the 
QR-I-CI calculations, MSD = 0.66%. Their agreement with the experimental data from 
[2611271128] is hne, too. Some larger deviations appear at the electron energies lower 
than 3 eV. Deviations from the theoretical data [29] are more noticeable. The reason for 
this discrepancy is that the BE, /-reduction procedure have not been applied in [29] . 

A completely different situation is observed for the 2s-3p excitation by electron 
impact. The cross sections of this process are presented in table [71 Our HF-I-CI results 
agree with the QR-I-CI results very well within all investigated electron energy range. 
Meanwhile, they are approximately by 34% smaller when compared to the data from 
[To] . It is important to notice that the 3p-2s radiative transition probability A values 
can change signihcantly, depending on the composition of the radial orbital basis and the 
amount of admixed conhgurations. Our calculated values converge to the data presented 
in table [5] only after the Cl expansion basis is made very large. At the same time, the 
electron-impact excitation cross section values have been rather stable and practically 
have not changed following the increase of the wavefunction basis, after all the most 
important admixed conhgurations have been included. 
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Table 6. The 2s-2p electron-impact excitation cross sections a (in 10 cm^) for Li 
atom at various incident electron energies e. 


e 

(eV) 

HF-hCI 

QR-fCl 

BEJ 

uni 

m 

Exp 

m 

m 

ccc 

m 

2.1 

18.459 

18.618 

18.374 


13.020 



2.3 

23.744 

23.884 

23.591 


17.595 



2.7 

30.236 

30.360 

30.017 


26.920 



3.5 

36.959 

37.065 

36.691 


35.717 



4.0 

39.258 

39.356 

38.978 


38.796 


36.08 

5.0 

41.762 

41.844 

41.474 


41.084 



5.4 

42.284 

42.360 

41.995 



49 


6.0 

42.738 

42.807 

42.451 




39.70 

6.6 

42.908 

42.971 

42.623 


41.172 



8.0 

42.611 

42.662 

42.333 




39.71 

10.0 

41.335 

41.373 

41.072 

38.00 


44 

38.33 

10.81 

40.689 

40.724 

40.433 


38.972 



15.0 

37.106 

37.126 

36.877 




34.47 

15.64 

36.570 

36.588 

36.346 


35.365 



20.0 

33.183 

33.194 

32.984 

33.10 


36 

31.36 

23.78 

30.655 

30.662 

30.473 


30.351 



25.0 

29.915 

29.921 

29.738 




28.87 

30.0 

27.221 

27.222 

27.061 




26.57 

38.6 

23.592 

23.591 

23.457 


23.691 



40.0 

23.096 

23.094 

22.963 




22.78 

50.0 

20.105 

20.101 

19.991 




20.09 

60.0 

17.839 

17.835 

17.739 

17.50 


28 


63.56 

17.160 

17.154 

17.064 


17.410 



70.0 

16.062 

16.057 

15.973 




16.35 

99.15 

12.532 

12.527 

12.465 


12.765 



100.0 

12.454 

12.449 

12.387 

12.40 



12.90 

149.4 

9.207 

9.202 

9.160 


9.387 



150.0 

9.179 

9.174 

9.131 

9.63 



9.60 

200.0 

7.333 

7.328 

7.296 

7.56 



7.66 

249.9 

6.139 

6.135 

6.110 


6.236 



300.0 

5.296 

5.292 

5.271 




5.56 

400.0 

4.183 

4.179 

4.164 




4.34 

400.5 

4.178 

4.175 

4.160 


4.230 



600.0 

2.982 

2.979 

2.970 




3.08 

601.4 

2.976 

2.973 

2.965 


3.010 



800.0 

2.337 

2.335 

2.330 




2.43 

802.3 

2.332 

2.330 

2.324 


2.348 



1000.0 

1.932 

1.930 

1.927 




2.01 

1404.2 

1.442 

1.441 

1.440 


1.447 



2000.0 

1.061 

1.060 

1.061 




1.08 
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Table 7. The 2s-3p electron-impact excitation cross sections a (in 10 cm^) for Li 
atom at various incident electron energies e. 


e 

(eV) 

HF-fCI 

QR-hCI 

BEJ 

uni 

4 

3.470 

3.486 

5.329 

5 

6.313 

6.317 

9.753 

6 

6.481 

6.481 

10.059 

7 

6.256 

6.254 

9.740 

8 

5.952 

5.950 

9.289 

9 

5.647 

5.645 

8.826 

10 

5.361 

5.359 

8.387 

15 

4.258 

4.256 

6.667 

20 

3.537 

3.535 

5.530 

30 

2.656 

2.655 

4.137 

40 

2.136 

2.135 

3.314 

60 

1.545 

1.544 

2.382 

80 

1.216 

1.215 

1.866 

100 

1.005 

1.004 

1.536 

120 

0.858 

0.858 

1.308 

140 

0.750 

0.749 

1.139 

160 

0.666 

0.666 

1.010 

180 

0.600 

0.600 

0.908 

200 

0.546 

0.546 

0.825 

250 

0.446 

0.446 

0.672 

300 

0.378 

0.378 

0.568 

400 

0.291 

0.291 

0.434 

600 

0.200 

0.200 

0.297 

800 

0.153 

0.153 

0.227 

1000 

0.125 

0.124 

0.184 

1500 

0.085 

0.085 

0.125 

2000 

0.065 

0.065 

0.095 

3000 

0.044 

0.044 

0.065 


We do not apply the /-reduction (IT^ for our data because the difference between 
our radiative transition probabilities and the data from [10] is only 6%, and the transition 
energies agree really well. Therefore, that reduction can not alter signihcantly our 
calculated electron-impact excitation cross section values. The overestimated cross 
section ctbej values in |T0] can be a consequence of the /-scaling procedure applied 
in that work. This is caused by the situation, where the El transition probability is 
proportional to the matrix element of the < > operator which is more important 

at longer distances. Meanwhile, when calculating the electron-impact excitation cross 
sections, one needs to apply the < fi{qr) > operator which has a maximum located 
closer to a nucleus. 

Unfortunately, we could not hnd any additional references on the experimental data 
or the calculation results of the 2s-3p excitation of Li atoms. 
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Figure 1. Comparison of the electron-impact excitation cross sections from the 

ground state to the 4p ^Pi /2 level. The dotted curve is the R-matrix (RM) damped 
cross section convoluted with a 30 eV Gaussian from [30], the solid curve represent our 
quasirelativistic (QR) calculation. 

f.f. Excitation of tungsten ions 

Previous comparisons for H, He, and Li demonstrate that the applied method is 
able to produce good-quality results for light neutral atoms. Described computer 
codes are mainly targeted for the production of electron scattering parameters for 
substantially more complex many-electron ions where it is very difficult to apply more 
accurate approximations and calculation methods developed for the electron-ion collision 
processes. In such a case, it is not a simple task even to determine correct and 
inclusive-enough Cl wavefunction expansion for the considered target system. The 
problems increase even more when one has to consider target-1-electron system containing 
signihcantly more conhgurations. This is the main reason why there are not so many 
data determined in R-matrix (RM) or convergent close-coupling (CCC) approach for 
heavy ions. Therefore application of rather simple approximation, such as plane-wave 
Born for the scattering process description becomes the most optimal, if not the only 
possible way so far to determine necessary data. 

There are just few data sources of collision data for complex many-electron tungsten 
ions. In [30] the results of fully relativistic, radiatively damped R-matrix calculations 
for and are presented. Authors in [30] were able to perform such kind of 

scattering calculations for these ions because of relatively small number of levels in the 
target expansion. We compare data from [30] to our calculation for the in Fig. [T] 

The threshold energy for this transition AE = 97.275 eV agrees very well with data 
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from [30] where Ai? = 97.233 eV. Similarly, the radiative transition probability in [30] 
A = 5.07 X 10^° s ^ is very close to onr determined valne oi A = 5.04 x 10^° s ^ (the 
deviation is jnst 0.6%). 

One can clearly see two different trends in that comparison. At the higher-energy 
end {E > 300 eV), the deviations are reasonably small, and they differ jnst by a 
few percent when the incident electron energy is above 800 eV. At the low incoming 
electron energies {E < 300 eV), the difference between onr data and R-matrix valnes 
is signihcant. The R-matrix calcnlation data are some 2.5 times larger than onr cross 
section valnes ai E = 125 eV, where cross sections reach maximnm valnes. Bnt that 
difference falls sharply when the incident electron energy increases. There are two main 
reasons for snch a behavionr. The R-matrix electron-impact excitation cross-section 
contains abnndance of resonances which reflect excitation throngh antoionizing levels of 
the ion. These resonances are located near excitation threshold, therefore they 

can make np 20% — 25% of the convolnted cross section valne, as one can see from figure 
4 in Another reason is that there is a strong interaction between incident free 

electron and the bound electrons of the ion. This interaction cannot be included 

in the plane-wave approximation adopted in the current work. Such an omission leads 
to smaller cross section values in our calculation. It is necessary to underline, that quite 
a large electron energy range is required when collision rates are determined. That leads 
to a rather small difference (especially at higher temperatures) in determined rate values 
applied in high-temperature plasma modelling. 

Even more complicated calculations are required when low ionization stages of 
tungsten ions are investigated. We have performed electron-impact excitation cross 
section calculations for several isoelectronic sequences such as Tm-, Yb-, Lu-, Hf- like 
ions including - W®’*' ions. Data for these ions are important in fusion plasma 
spectra modelling. We have determined both the spectroscopic parameters (energy 
levels, transition wavelengths, radiative transition probabilities) and the electron-impact 
excitation parameters (collision strengths, cross sections, collision rates) for several ions 
of these sequnces. Our data are incorporated in Atomic Data and Analysis System 
(ADAS) [5] as the basic parameters and are utilized to determine the derived parameters 
for various plasma conditions. 

To demonstrate the abilities of our developed computer codes, we present several 
plots of the electron-impact excitation cross sections for the tungsten ions ranging from 
W®’*' to W^"*". Our performed calculation involves a large number of excitations involving 
levels of these ions. All possible excitation transitions involving several hundreds of levels 
are considered. Plots given in the present paper are given for demonstration purpose. 

Figures [2] and [3] demonstrate cross sections for the W®’*' ion excitation from the 
ground and from the first excited level. For other three ions, namely W^+, W^^ and 
W^"*", we present only the excitation parameters from the ground level in Figs. 01 E] E] 
although much more data are produced. For example, we consider the lowest 293 levels 
in W^"*" and transitions among them. 


Cross section (10’'® cm^) Cross section (10’'® cm^) 
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Figure 2. Electron-impact excitation cross sections (in 10“^® cm^) from the ground 
state 5p®5d ^ 03/2 of the W®"*" ion to the excited states 5p®5d ^ 05/2 (solid curve), 
5p®6s ^Si /2 (dashed curve), 5p®6p ^Pi /2 (dotted curve) and 5p®6p ^P 3/2 (dash-dotted 
curve). Incident electron energies are given in transition threshold units. X = E/AE. 



Energy (threshold units) 


Figure 3. Electron-impact excitation cross sections (in 10“^® cm^) from the first 
excited state 5p®5d ^ 05/2 of the ion to the excited states 5p®6s ^Si /2 (solid 

curve), 5p®6p ^Pi /2 (dashed curve) and 5p®6p ^P 3/2 (dotted curve). Incident electron 
energies are given in transition threshold units. 











Cross section (10'^® cm^) Cross section (10'^® cm^) 
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Figure 4. Electron-impact excitation cross sections (in 10“^® cm^) from the ground 
state 5p®5d^ ^F 2 of the ion to the excited states 5p®5d^ (solid curve), 

5p®5d^ ^Pq (dashed curve), 5p®5d6s ^Di dotted curve) and 5p®5d6p ^F 2 (dash-dotted 
curve). Incident electron energies are given in transition threshold units. 



Energy (threshold units) 


Figure 5. Electron-impact excitation cross sections (in 10 cm^) from the ground 
state 5p®5d^ "^F3/2 of the ion to the excited states 5p®5d^ "^F 5/2 (solid curve), 
5p®5d^6s '*F 3/2 (dashed curve) and 5p®5d^6p ^F 3/2 (dotted curve). Incident electron 
energies are given in transition threshold units. 
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Figure 6. Electron-impact excitation cross sections (in 10“^® cm^) from the ground 
state 5p®5d^ ®Do of the ion to the excited states 5p®5d^ ®Di (solid curve), 

5p®5d^6s ®F 2 (dashed curve) and 5p®5d^6p ^Di (dotted curve). Incident electron 
energies are given in transition threshold units. 


5. Summary and conclusions 

The performed calculations of the electron-impact excitation data and their analysis 
have demonstrated that the computational methods implemented in our newly- 
developed computer codes are suitable to produce the electron-impact excitation cross 
sections and related parameters, such as the electron-impact collision strengths or 
collision rates. One can expect that such the data for other ions will be reliable enough. 

Based on good agreement of the radiative transition probabilities with existing 
data, we expect that our determined electron-impact excitation cross sections for the 
2s-3p transition in Li atom are accurate enough. The comparison of cross sections 
for leads to conclusion that the calculated parameters even for highly-charged 

ions can be applied in high-temperature fusion plasma spectra modelling. Signihcant 
deviations from earlier published results can be explained by the fact that radiative 
transition partameters are influenced by inclusion of correlation corrections differently 
compared to excitation cross sections. Therefore a simple normalization described by 
Eq. flT2l) in [lO] is not well-grounded. 

We are planning to incorporate results of such calculations (together with other 
spectroscopic parameters) into the newly-developed database ADAMANT when the 
electron-impact excitation processes in various many-electron atoms and ions with open 
S-, p-, d- and f- shells are considered. 
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